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A  Galerkin  Method  on  Nonlinear  Subsets  and 


Its  Application  to  a  Singular  Perturbation  Problem 


ABSTRACT 


In  the  Ritr-Galerkin  method  the  linear  subspace  of  the  trial  solutions  is  extended  to  a  closed 
subset.  As  an  example,  a  class  of  so-called  subiinear  approximation  and  interpolation  is  developed. 

Some  results,  such  as  orthogonalization  and  minimum  property  of  the  error  function,  are  obtained. 

*» 

A  second  order  scheme  has  been  developed  for  solving  the  linear  singular  perturbation  elliptic 
problem . 

-  <  u”  +  p(x)  u’  +  q(x)  u  -  f(x),  u(0)  «u(l)  -  0. 

Error  estimates  are  given  for  a  uniform  mesh  size  h: 

IK1  •  «lli  i  c,  k1-5.  n.™,11  -  «,i,||00  i  c^.  hw,  O—o.i) 
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if  h  < -  c,  where  the  constants  C-  and  C  .  (i— 0,1)  all  are  uniformly  bounded  for  small  c. 

IIPlI 

00 

For  the  same  accuracy,  the  present  nonlinear  scheme  is  one  order  of  magnitude  more  than  the 
usual  method  used  in  the  piecewise  linear  subspace.  Numerical  results  for  linear  and  semi-linear 


singular  perturbation  problems  are  included. 
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1.  Introduction 

The  development  of  finite  element  methods  has  been  successful  in  various  fields.  From  a 
mathematical  point  of  view,  the  method  is  one  of  extensions  of  Rayleigh-Ritr-Galerkin  technique, 
([1],  [15],  [16],  [17]).  Usual  finite  element  schemes,  choosing  piecewise  polynomials  as  trial  functions, 
are  very  efficient  when  there  are  no  steep  gradients  in  the  true  solution.  Otherwise,  poor  results 
might  occur.  In  order  to  get  accurate  numerical  data,  one  may  use  adaptive  mesh  technique(e.g.  [8]) 
or  a  higher  precision  scheme  such  as  h-version  and  p-version  respectively  [3].  Beyond  usual 
polynomials,  rational  elements(e.g.  [24])  and  exponential  elements  [9]  have  been  introduced  to  enrich 
the  trial  subspace  to  reduce  number  of  parameters  for  a  given  precision.  One  thing  in  common 
among  these  techniques  is  that  they  are  all  reduced  to  a  discrete  linear  system  if  the  original 
differential  equation  is  linear. 

Nevertheless,  our  approach  is  quite  different.  To  find  a  better  discrete  approximation  of  weak 
solutions  with  steep  gradients,  we  try  to  relax  the  limitation  of  replacing  the  continuous  variational 
problem  only  by  a  sequence  of  finite-dimensional  subspaces.  Hence,  in  this  paper,  we  present  an 
extension  of  the  finite  element  method  from  subspace  to  more  general  subsets  and  adopt  the  method 
to  solve  singular  perturbation  problems  (including  linear  and  semi-linear)  in  one  dimension.  For 
linear  problems,  our  aim  is  to  solve  a  small  semi-linear  system  instead  of  a  large  linear  system  which 
arises  by  using  the  usual  trial  subspace  of  piecewise  polynomials  for  a  given  precision. 

From  a  practical  point  of  view,  there  are,  at  least,  two  questions  which  need  to  be  answered 
now.  First,  how  to  find  a  good  non-linear  approximation  of  a  non-linear  functional  space  which  can 
be  devised  especially  for  singularity  problems.  Secondly,  how  to  solve  the  resulting  discrete  non¬ 
linear  system  efficiently.  This  non-linear  approximation  should  include  conventional  piecewise 
polynomial  and  it  is  expected  to  be  not  too  far,  in  some  sense,  from  linear  approximation  in  order  to 


meet  the  theoretical  demand  (such  as  convergence  and  to  keep  some  behaviors  of  the  true  solution) 
and  to  satisfy  the  practical  aim. 

The  approximation  used  in  this  paper  is  called  pitetwiat  mapping-polynomial  or  apline 
mapping-polynomial.  It  means  that  the  approximation  is  of  piecewise,  and  in  each  subinterval  a 
local  one  to  one  mapping  is  applied  first,  then  a  polynomial  approximation  is  used  in  the  mapping 
plane.  The  final  approximation  is  obtained  by  using  the  inverse  mapping,  and  the  whole 
approximation  function  has  some  orders  of  smoothness  according  to  various  requirements.  In 
particular,  it  reduces  the  usual  polynomial  or  polynomial-spline  approximation  if  the  mapping  is 
always  equal  to  the  identity  mapping. 

A  large  amount  of  attention  has  recently  been  focused  on  the  difficult  singular  perturbation 
boundary  value  problems.  These  problems  arise  from  some  different  fields,  for  examples,  boundary 
layer  or  convective-diffusion  type  flows  in  fluid  dynamics.  Conventional  methods  applied  to  such 
problem  result  in  unrealistic  oscillation  and  poor  approximation  unless  the  mesh  length  h  is 
excessively  small.  Some  effects  have  been  done  by  various  authors  using  local  higher  order 
polynomial  approximation  with  some  parameter,  called  ’Upwinding'  methods,  to  match  the  true 
solution  better  at  the  nodes.  The  method  has  been  discussed  by  Christie  and  Mitchell  [0], 
Barrett, Morton  [4],  Heinrich  and  el.  (12],  Babyska  [2],  etc.  An  "Exponentially  fitted  method" 
developed  by  de  Groen  and  Hemker  ]9j.  is  to  add  a  piecewise  exponential  term  to  enrich  the 
subspace  of  piecewise  polynomial. 

In  section  2  and  3,  we  generalise  respectively  the  usual  Ritz  and  the  Gaierkin  method  from 
linear  trial  subspaces  to  subsets,  and  derive  some  results  such  as  orthogonalization  and  error 
estimations.  A  brief  discussion  about  ’sub-linear'  operator  and  its  approximation  is  given  in  section 


4 


4.  In  section  5,  the  semi-linear  finite  element  technique  is  studied  by  solving  singular  perturbation 
problems  in  one-dimension:  -€u”+pu’+qu“f,  u(0)=»u(l)=»0.  The  results  show  an  improvement  over 
one  more  order  precision  than  the  corresponding  scheme  of  using  piecewise  linear  subspace  and  that 
the  constraint  of  mesh  size  h  is  relaxed  from  0(<2)  to  0(e).  A  linear  and  semi-linear  test  singular 
perturbation  problems  are  given  in  section  6.  Computational  result  agree  with  the  above  theoretical 
analysis. 

Some  research  results  on  the  same  topic  in  two-dimensions  will  be  reported  separately  [22]. 
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2.  A  Rits  method  on  subsets  for  self-adjoint  equations 

Consider  a  self-adjoint  elliptic  linear  differential  equation 

Lu  -  f  (l) 

a(u,v)  —  (Lu,v)  is  a  positive  quadratic  form  in  a  real  Hilbert  space  H  with  an  inner  product  (*,*)  and 
a  norm  ||*||: 

c2  Hull2  a(u,u)  C,  ||u||2,  for  all  ucH  (2) 

where  and  C2  are  positive  constants.  It  implies  that 

|a(u,v)|  £  C,  ||u||  ||r||,  for  all  u,v  t  H.  (3) 

u  is  defined  as  a  weak  solution  of  (1)  if  it  satisfies 

a(u,v)  — >  (f,v)  for  all  v  <  H.  (4) 

It  is  well  known  that  u  is  a  weak  solution  of  (1)  if  and  only  if  u  is  the  unique  minimum  solution  of  a 
quadratic  functional  I,  i.e. 

I(u)  —  inf  I(v)  —  inf  {a(v,v)  -  2(f,v)}  ~  (5) 

«H  »«H  '  ' 

As  a  well-known  discrezation,  H  in  the  variational  problem  (5)  is  replaced  by  a  sequence  of 
finite-dimensional  subspaces  Vh  contained  in  H: 

I(uh)  -  mf  I(v) 

TfV* 

which  is  equivalent  to  the  following  weak  solution 

a(uh,  vh)  —  (f,  vh),  for  all  A  c  Vb.  (g) 

Now  we  replace  H  in  (5)  by  a  sequence  of  closed  subsets  Sh  with  the  same  finite-dimensional 
parameters.  Let  T  be  an  one-to-one  continuous  mapping  from  an  open  convex  set  of  V*1  onto  S*1: 
TV,h  -  Sh. 

Definition  1:  [25]  The  mapping  T:  Vjh  — >  S*1  is  differentiable  in  the  open  convex  set  vA 

if  for  each  v«Vj  there  is  a  Jacobian  matrix  T’(v)  such  that 
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Hm  ||- 

tt->0 


v+an  -  T  v 


T'(v)q||  «*  0,  for  each  q  t  Vj 


In  particDlar,  T’  —  T  if  T  is  a  linear  mapping. 

Consider  a  restricted  variational  problem  on  the  closed  subset  S*1: 

I(u  )  —  inf  I(v)  (8) 

8  ™s‘ 

Since  Sh  is  closed,  so  there  exists  a  solution  of  (8)  in  S1*.  If  ug  minimizes  I  over  Sh,  ug*Tw,  then  for 
any  a>0  and  q  <  Vjh 

I(u)£.I(T(w+aq)).  (0) 


T(w+aq)  «■  Tw  +  aTq  +  ic(a) 

where  T  is  positive-homogeneous  and  — 

#c(a)  ■»  T(w+aq)  -  Tw  -  aTq 
The  right  side  of  (9)  is 

I(ug)  +  2a|a(ug,Tq)  -  (f,Tq)]  +  2[a(ug,ie(o))  -  (f,ic(or))] 

+  a2a(Tq,Tq)  +  2aa(Tq,fc(a))  +  a(ie(a),»e(a))  —  1(a) 

As  a  function  of  the  parameter  a,  the  fact  that  u.  minimizes  I  over  S  requires  lim  I’(a)  ■»  0. 

8  o->0 

Observing  that 

«(o)-o,  K’(o)-(r(r1ug)-T)q, 


0  “  I’(«)la_>0  -  2  ( «(Vr*)  -  ^ >  -  &  K'(° ) )} 

hence,  it  yields 

a(ug,T’(T1ug)q)  -  (f^r1^)  for  all  q  <  V,h. 

Therefore 

Theorem  2:  If  (i)  V*  is  a  subspace  of  H.  (it)  S*  is  a  close  subset  of  H.  (Hi)  T  is  an  one- 


to-one  positive  homogeneous  and  differentiable  mapping  from  an  open  convex  set  V}  of  V* 
onto  Sk:  TVj  «  5*.  Then  (iv)  There  exists  a  solution  uf  of  (8)  and  (10)  holds. 

The  above  Theorem  shows  that  the  nonlinear  system  (10)  has  at  least  a  solution  which 
minimises  the  variational  problem  (8).  Usually,  it  does  not  mean  they  are  equivalent  each  other. 
Because  there  are  no  guarantee  of  unique  solution  in  general  case.  However,  we  have  the  following 
conclusion: 

Theorem  3:  If  Vj  contains  uh  defined  in  (6),  then  for  the  mapping  T  which  is  sufficient 
close  to  a  linear  mapping,  i.e.,  ||T  -  T’  ||  is  sufficient  small,  the  nonlinear  system  (10)  has 
unique  solution  which  minimizes  the  variational  problem  (8). 

Proof:  In  fact,  (10)  can  be  rewritten  as 

a(u,  v*1) —  (f,  vh)  +  Q^v*1) 

where 

QKv**)  -  a(u,  [T-T’(r1u)]vh)  +  (f,  [T-T’(r1u)]vh)  )  — 

Since  there  is  unique  solution  in  (6),  hence,  the  above  equations  system  also  has  unique 
solution  if  ||T  -  T’||  is  sufficient  small. 


Now  we  suppose  that  the  generalised  coordinates  (real  parameters)  of  the  subset  Sb  are  q^...^, 
then  the  first  variational  equations  of  I(w)  in  Sh  must  be  vanished 

\  glj.  ™  °»  for  '  ™  (n) 

and  the  determinant  of  the  second  variational  matrix  at  the  point  of  the  solution  is  positive 


(12) 


Let  {Bj}  be  a  basis,  then  for  each  w<Sb 

w  ■>*  T*w  +  w*,  ^  —  B.  +  ,  where  T*w  “  E  qjBj,  w*  — «  w  -  T*w. 


Substituting  the  above  formulas  into  (11)  yields 
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5 


p 

>  ■ 

z ; 


r' 


4 


S 


4 


4 


4 


Ea(B.,Bj)qj-(f,Bi)  +  G.(q) 
where  G.  —  (f,|~)  -  *("*  J““)  *  £ 


(13) 


Hence,  the  equations  of  the  weak  solution  in  subsets  are  different  from  ones  in  subspaces  only  by 
the  last  extra  term  which  tends  rero  when  the  subset  Sk  tends  a  subspace.  Also,  the  system  (11)  can 
be  written  as 


a(w,B.)  -  (f.Bj)  +  G*j(q),  where  G*j  -  (f.^1)  -  a(w,§2L). 

Hence,  for  each  vfVk,  ignoring  the  extra  terms,  we  get  an  approximate  equations 
a(ug,v)  «=  (f,v)  for  all  veVh 


(14) 


(15) 


Because  w  in  (14)  corresponds  to  the  unique  solution  of  the  variational  problem  (8)  for  the 
positive  quadratic  form  a(u,u)  restricted  in  the  subset  Sk,  being  the  continuity  of  solutions  with  the 
system,  there  also  exists  a  solution  ug  of  the  system  (15)  in  Sh,  if  the  distance  between  Vk  and  Sk  is 
sufficient  small.  Geometrically,  it  is  obvious.  In  fact,  from  (11)  and  (12),  it  means  that,  as  a 
hypersurface  in  the  n  dimension  of  (q|,...q„),  *  *  dl/dqi  is  separated  by  a  hyperplane  z  «■*  0  and  they 
have  only  one  intersection  point.  Moving  this  hypersurface,  there  still  exists  a  unique  intersection 
point  if  the  moving  distance  is  sufficient  small. 

For  practical  aim  there  is  another  approximation  versions  of  (14):  Find  ugtSh  such  that 

a(ug,  ug  -  v*1)  -« (f,  ug  -  v1*),  for  all  v*  «  Vh.  (16) 

Suppose  ug  is  the  unique  solution  of  (16).  From  (4),  for  any  v*1  in  Vh,  a(u,ug  -  v*1)  -«  (f,ug  -  vh), 
subtracting  (16)  leads  to  a(u  •  ug,ug  -  v1*)  0.  Hence 

a(u  -  vh,u  -  vh)  =■  a(u  -  ug,u  -  ug)  +  afv*1  -  ug,vh  -  ug). 

Using  (2),  furthermore,  for  any  v11  in  Vh, 
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Co  ||u  -  ug||2  a(u  -  ug,u  -  ug)  a(u  -  vh,u-vh)  C,  ||u  -  vh||2. 

There  are  similar  formulas  for  the  case  of  (15).  Thus,  we  have  proved  the  following  fundamental 

theorem  of  the  Ritz  method  on  subsets  which  is  an  extension  of  the  Theorem  1.1  in  [2]  for  subspaces. 
Theorem  4:  Suppose  ug  is  the  unique  solution  of  (16)  or  (15)  in  a  closed  subset  5*,  then 
it  satisfies  the  following  properties: 


(a )  Minimum  property 

a(u  -  u  ,u  -  u  )  —  inf  a(u  -  Au  -  v11), 
s  s  Tk<vk 

or 

a(u  -  us,u  -  ug)  —  inf  a(u  -  ug  -  vh,u  -  ug  -  v*1), 
and 

lh-  us|(  ^  C  jnf  ||u  -  vh|| 
or 

ll«  -  «gll  C  jnf  ||u  -  ug  -  ^11 
where  C  is  a  constant. 


(17) 

(18) 


(19) 

(20) 


(b)  Orthogonalization 

a(u  -  ug,ug  -  v*1)  0,  for  all  v*1  in  V*1.  (21) 

or 

a(u  -  ug,vh)  **  0,  for  all  vh  in  Vh.  (22) 


In  practical  view,  as  a  system  for  the  weak  solution,  (15)  is  more  attractive  than  (16).  And  the 


difference  between  them  could  be  small  if  the  subset  is  'not  far’  from  a  subspace  in  some  sense. 


3.  A  Galerkln  Method  on  a  closed  nonlinear  subset 


The  analysis  in  section  2  can  be  extended  from  the  Rite  to  the  Galerkin  method.  Suppose  that 
the  operator  L  in  (1)  is  not  self-adjoint  in  which  derivatives  of  odd  order  spoil  the  self-adjointness  of 
an  elliptic  equation  and  the  associated  quadratic  functional  I(v)  defined  in  (5)  is  not  positive  definite. 
The  problem  now  is  to  And  a  stational  point  rather  than  a  minimum  of  I(v).  There  are  some  results 


on  the  existence  of  the  weak  solution  (4),  e.g.  Babuska  and  Asia  [1],  Strang  and  Fix  [17].  Let  us 
quote  a  few  results  of  Galerkin  Method  first. 

Theorem  S:  Let  Hj  and  Hg  be  two  real  Hilbert  spaces  with  inner  products  and 

(*,  *)fj ,  respectively,  (f,v)  be  a  continuous  linear  functional  on  Hg  and  a(u,v)  a  bilinear 
form  with  three  inequalities 


(i)  |a(u,v)|  Cj||u||H  ||v||j^  for  all  otH,  and  ▼  t  H^,  where  C,  <  oo. 
|a(u,v)| 

(ii)  inf  sup  .  C2  >  0  . 

««.  «h,  Hulln  lMii^ 

(iii)  sup  |a(u,v)|  >  0,  v  yt  0 


Then  there  exists  one  and  only  one  weak  solution  ug  of  the  functional  equation  Lu 
such  that 


f 


a(«0>v)  —  (f»v)  for  all  v  <  Hj 

and 


(23) 


IKIIh,  ^  C21  llflltf, 

A  proof  of  this  result  can  be  found  in  [1],  theorem  5.2.1.  Galerkin’s  method  is  the  natural 

discretization  of  the  weak  form.  In  general  it  involves  two  families  of  functions _ a  subspace  Sh 

of  the  solution  space  (or  trial  space)  Hj  and  a  subspace  V*1  of  the  test  space  Hj.  Then  the  Galerkin 
solution  uh  is  the  element  of  Sh  which  satisfies 

a(uh,  vh)  -  (f,vh)  for  all  v*  <  Vh  (24) 

Since  both  S*1  and  Vh  are  linear  subspaces,  if  (sj)  is  a  basis  for  Sh  and  { )  is  a  basis  for  Vh,  the 


(25) 


solution  un  «  Eq^Sj  satisfies  a  linear  system 
Aq«d 

where 

A  ««  (  afsj.vp  ),  d  —  (f.Vj)  (20) 

If  A*1  exists,  there  is  a  unique  solution  ub  of  (24).  Also,  there  are  some  error  estimations  of  the 
Galerkin  method,  say,  see  Strang,  Fix  [17]  and  Aziz  [lj.  However,  if  there  is  an  odd-derivative  term 
of  the  bilinear  form  with  significant  size,  the  Galerkin  method  is  usually  unsatisfactory.  The 
essential  reason  is  that  the  approximation  in  linear  functional  space  is  not  good  enough  in  this 
singular  case.  Probably,  that  is  one  way  to  overcome  the  difficulty  is  to  extend  the  trial  solution 
space  to  a  nonlinear  subset. 

Now,  suppose  that  Sh  which  is  a  closed  subset  of  H}  has  the  same  number  of  freedoms  with  Vh 
and  that  there  exists  a  element  uh  <  Sh  such  that  (24)  still  holds  true.  Being  (23),  subtraction  yields 
the  following  Lemma. 

Lemma  0:  For  any  subset  of  Hv  if  there  exists  an  element  u*  t  Sk  which  satisfies  the 
relation  (24),  then  with  respect  to  the  energy  inner  product,  u*  is  the  projection  of  u  onto 
5*,  or,  the  error  u  -  uk  is  orthogonal  to  V* 

a(u  -  uh,  v11  )  — «  0  for  all  v1*  e  Vh  (27) 

Let  the  notation  Uj  denote  an  interpolation  of  any  u  c  in  the  subspace  Vb.  Since  for  any  Uj  c  Sb, 
a(u  -  uh,  u  -  uh)  —  a(u  -  uh,  u  -  Uj)  +  a(u  -  uh,  Uj-uh) 
being  (27), 

a(u  -  uh,  uruh)  —  a(u  -  uh,  (uj-uhHuruh)j) 
or 


a(u  -  uh,  Uj-uh)  ■■  a(u  -  u\(u  -  uh)-(u  -  uh)j)  -  a(u  -  uh,  (u  -  UjM«  *  «j)j) 


So,  from  (he  inequalities  of  Theorem  5 


C2||u  -  uh||2Hj  ^  C,||u  -  uh||Hi{||u  -  Ujll^  +  ||(uj-®hM®j-oh)1||H2} 

Therefore,  we  proved  following  error  estimations. 

Theorem  7  s  Suppose  the  conditions  in  Theorem  5  hold  as  well  as  (24),  then  on  the  closed 
nonlinear  set  Sh  the  approximation  solution  u*  of  (27),  if  it  exists,  has  following 
estimates 


ll«  -  ah||Hi  Hu  ‘  uh  ‘  wllj^ 

II®  -  ®hllHi  (II®  -  “jIIhj  +  IKuruhHvuhVy 


(28) 


or 


n®  -  «hiiHl  &.-7T  ( ii(®*®h)  ■  («-«h)iiiH2 + + ik®-®j)  •  (29) 


BjcS1 

dr 

Corollary  1.  If  the  subset  Sh  coincides  with  the  subspace  Vb,  then 
C, 


II®  -  ®h||H  ^  ~  »®{  II®  -  ® jIIil  (30) 

J  C2  «J<5‘  *  "a 

(30)  is  just  the  result  of  the  usual  Galerkin  method.  Hence,  (24)  above  is  just  a  generalization  of  the 
Galerkin  method. 


Corolla  7  2.  Let  uJh  be  an  interpolation  of  u  on  the  subset  Sh,  then 

II®  •  ®hllHl  ^  ^(11® *  «jhllH2  +  IK*  *  ahM»  *  AHh,  +  +  IKU  *  uJhH®  -  “jhVlHjI  (31) 

The  bounds  (29)  •  (31)  will  play  a  central  role  in  error  analysis  .  It  is  clear  that  the  subset  Sn  should 
be  so  chosen  as  it  can  tends  a  denumerable  dense  set,  as  h  tends  zero,  in  the  true  solution  space  H,, 
as  well  as  Vh  in  H2.  In  this  case  the  limiting  behaviors  of  the  error  in  energy  norm  as  h~>0  depends 
mainly  on  the  approximation  ability  of  the  subset  Sk. 

Now  we  turn  to  discuss  existence  and  uniqueness  of  solution  of  (24)  briefly.  It  was  considered  in 
section  2  for  self-adjoint  a(u,v).  When  a(u,v)  is  not  self-adjoint,  in  terms  of  variational  principles, 
the  weak  solution  (24)  is  equivalent  to  find  the  stationary  point  for  the  bilinear  a(u,v)  on  Sk  x  Vh 
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where  StHj  and  VcH„  are  a  closed  approximation  subset  and  a  subspace  with  same  finite  parameters 
respectively.  Because  the  existence  of  the  stationary  point  in  the  whole  space  Hj  is  assured  by 
theorem  5,  hence,  from  geometric  intuition,  there  is  a  stationary  point  in  the  sense  of  (24)  for 
sufficient  small  h,  at  least.  Besides,  if  there  is  unique  stationary  point  of  (24)  when  the  subspace  Sk 
coincides  with  the  subspace  Vb,  then,  the  stationary  point  still  exists  if  the  subset  Sh  is  ’very  close’  to 
the  subspace  Vh.  In  general,  we  have 

Theorem  8:  Suppose  there  exists  a  svbspace  SLh  with  a  basis  («y)  in  which  the  linear 
system  (24)  has  unique  solution  and  T  is  a  map  from  the  subset  S1*  to  the  subspaee  SLh 
such  that  for  a  basis  (ik)  of  the  test  subspaee  V* 

p  (A-*J(G) )  <  1  (32) 

where  the  notation  p  denotes  the  spectral  radius  of  a  matrix,  A  defined  in  (26  ),  and  J(G) 
is  the  Jacobi  Matrix  o f  the  vector  G  defined 

G  **  (  a(uh  -  Tuh),  v.  ), 

then  the  nonlinear  system  (27)  exists  unique  solution. 

Proof;  Let  Tuh  —  E  q-.v,  since  a(uh,v.)  —  a(Tuh,v.)  -  a(uh  -  Tuh,Vj)  ,  from  (24),  (27) 
becomes  E  a(sj,Vj)  «■=  (f,Vj)  +  Gj,  In  matrix  form  it  can  be  written  as 

Aq*d  +  G(q).  (33) 

Using  the  following  ’simple’  iterative  procedure 

A  q<°>  -  d, 

A  qM  -  d  +  G(q<k'1))  (34) 

which  is  a  contraction  mapping  if  the  condition  (32  )  is  satisfied.  Q.E.D. 

Remark:  (33)  is  very  useful  not  only  for  proving  existence  of  the  solution,  but  also  for 
computing. 

For  practical  view,  hence,  the  first  problem  for  using  the  generalized  Galerkin  method  is  to 
construct  an  adequate  nonlinear  approximation  subset  as  Sh  above. 
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4.  ’Sub-linear’  approximation  and  interpolation 

Let  T(u)  be  a  real  operator  of  u,  where  u(x)  be  a  real  function  defined  a  given  vector  space  X 

and  belong  to  a  space  S,  T(u)  belong  to  S,  too. 

Definition  9:  An  operator  T(u)  is  called  positive  on  the  set  X,  if 

T(u)  >  0  for  all  u(x)  >  0  and  x  t  X  (35) 

Definition  10:  An  operator  T(u)  is  called  sublinear  on  the  set  X  if  it  satisfies  two 
following  conditions 

(i)  Positive-homogeneous 

T(au)  =  aT(u)  for  all  a  >  0  in  R  and  u  <  U,  x  t  X  (36) 

(ii)  Subadditive 

T(u+v)  ^  T(u)  +  T(v)  for  all  u,v  <  U  and  x  <  X 
or 

T(u+v)  ^  T(u)  +  T(v)  for  all  u,v  e  U  and  x  c  X.  (37) 

Consider  interpolation  and  approximation  using  sublinear  piecewise  positive  operator.  For  simplicity, 
let  the  set  X  ■■  (0,1],  and  a  partition  A  be  given 

A:  0  -  x0  <  x,  <...  <  xN  -  1,  hj  -  Xj  -  x..,  (38) 

Particularly,  the  linear  positive  operator,  defined  by  Korovkin[l3]  is  sublinear  positive. 

When 

Bj(x)  ^0,  E  Bj(x)  «  1,  for  all  x  in  [0,1]  (39) 

then 

T(u)  —■  E  u(x.)  Bj(x)  (40) 

is  positive.  As  an  example,  can  be  chosen  as  B-spline.  Similar,  if 

B(t,x)  2. 0,  for  all  tpc  in  (0,1],  and  J0l  B(t,x)  dt  1,  for  all  x  in  (0,1].  (41) 

then 


(42) 


T(u)  =  /0>u(t)  B(t,x)  dt 
is  positive  too. 

Lemma  11:  If  u^  Vy  >  0,  p  >  1,  Bjfx)  is  defined  by  (89),  then 

{E(uj+vj)PBj(x)}1/P^{E(uj)PBj(x)}1/P+{E(vj)PBj(x)},/P.  (43) 

The  inequality  direction  will  be  opposite  ifp<l(py^0).  In  the  limit  ease  of  p  «=  0, 

(48)  becomes 

ntUj+VjfjW  +  I1(t.*|M.  (44) 

In  each  case  the  equality  holds  true  if  and  only  if  the  two  sequences  (  u  )  and  (  v  )  are 
proportional. 

Id  fact,  the  above  inequality  is  just  the  triangular  inequality  tor  the  lp  space  with  weight.  It  can  be 
easily  proved  using  a  classical  inequality,  e.g.,  {5}.  Hence,  the  operator 

T(u;p)-  {E(uj)PBj(x)}1/*>  (45) 

is  sub-linear  positive. 

For  instance,  if  we  take  the  basic  functions  (Bj(x)}  as  B-spline  and  u.  as  an  average  of  u(x)  on 
some  nodes  near  Xj,  then  (45)  becomes  a  sublinear  positive  approximation  operator  of  u(x)  on  [0,1],  it 
wll  be  a  generalization  of  the  well-known  Schoenberg  approximation. 

Consider  a  kind  of  piecewise  interpolations  using  the  above  semi-linear  positive  operator.  For 

Xj.,  x  <L*y  let  t  —  (  x  -  Xj.,  )/  h.,  ^  —  o(t)|^0,  u,  —  u(t)|t_,,  p  —  Pj,  and 

T(u;  p)  —  {  uf  (1-t)  +  u,p  t  }1/p(pj*0)  0  <  t  <  1  —  (46) 

T(u;  0)  —  u^1*^  u,1  (p  —  0)  0  <  t  <  1  (47) 

Obviously  T(u;  p)  is  piesewise  sublinear  and  positive. 

Theorem  12:  The  interpolator  operator  T(u;  p)  is  piecewise  sublinear  and  positive,  and 
if  u  iC^lOyl),  for  {  <  z  <  Xj  then  for  u(z)  >  0  there  is  a  remainder  expression 

u(x)-T(u;p)  —  5(x-xj.1Xx.-xXo”*(l-Pj)«’2/»)lx.  ^+0(hs). 


id 


i 


i 


i 


> 


l 


i 


(  Xj.,  <  {  <  X.  ) 

(48) 

Furthermore 

Maxju(x)-T(u;  p)(  ^-Max|u”-(l-Pj)u'2/u|  +  0(h3). 

(49) 

and 

Max|u’(xXT’(u;  p)|  ^^Max|u”-(l-Pj)u’2/u|  +  0(h2). 

(50) 

Proof:  Since  T(u)“U  for  t»0  and  t— <1,  using  a  well-known  technique  of  error  estimates  in 
Lagrange  interpolation  leads  to  (48)  directly,  so 


u(x)  -  T(u;  p)  —  ^x-Xj.jXxj-xX  n”-T’’)|xi_  f  ,  (xj.,<e<xj). 
from  the  Taylor  expansion  T”  =*  (l-pj)u’2/u  +  0(h2).  Hence,  (49)  and  (50)  follow. 

In  particular,  if  p  — «  1  everywhere,  (46)  becomes  piecewise  linear  interpolation.  Hence,  in  general 
case,  now  the  piecewise  linear  operator  operates  on  up  instead  of  on  u  itself.  In  another  words,  (46)  is 
a  generalized  means  replacing  the  arithmetic  means  with  weights  for  the  linear  case  (Jiachang  Sun 


(19]).  Furthermore,  if  we  choose  the  piecewise  constant  Pj  such  that 


the  resulting  interpolation  (46)  will  have  one  more  order  of  precision. 

Theorem  13:  If  v  t  C^fO,l],  then  for  the  pieeewiee  interpolation  (46)  with  (51) 

Max|uW  -  T(i)(u)|  £  Cihs’iMax|W(u)|  +  Oft4-1).  i— 0,1,2. 

||u  •  T(u)||l,  Ch3Max(W(u)|  +  0(h4) 
where 


ol/2  ill  n’2  u 

co“  ci“  i?  c2  -  7- c  -  w(u>  “ 

Proof:  Set  p«—p,.  Applying  the  Taylor  expansion  yields 

UU”  X.+Xfl  u  . 

p-1 - ^+(-L_^..x)(-^”  +  0(h2), 

u  2  2  u  ^ 

(uP-^’V  -  uP-2u’2(^l.  x)  (1)”  +  0(h2), 

2  u’ 

(uP-V)”  -  -uP'2u,2(^-)1’  +  0(h). 

u’ 

Using  these  results,  a  straightforward  computation  leads  to 


(51) 

(52) 

(53) 
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«,p  “op  °p  +  ^x^xXx-xoKop■1®,), 

+  ^Xj-xKx-XoKxj^x+x^uP-^’)”  +  0(h4) 

-  »p  +  iWxXx-x0X^  -  +  0(h4). 

6  2  u 

and 

j  pt  2 

-(up-,u,)”{(xI-x)2  -  (x,-xXx-x0)  +  (x0-x)2>  +  0(hs) 

6 

-  uP'1u’{l+1ju(^”[(xrx)2+(x0.x)2-4(x1.xXx-x0)]}  +  O (h*) 

Hence 

x^  x.-x  , , 

T(u;p)  —  (  „0P  _L_)1/P 

l  ^  |  j.  ^ 

-  n  {  1  +  ^xrxXx-x0X— -  xK2u’2(-^”}  +  0(h4). 

Therefore 

T(u;p)  -  «(x)  -  ^Xj-xXx^X^-  -  x)^y’’  +  0(h4). 

Since  the  function 

h'3|J(xrxXx-x0X^2.  .x)|-|^l-t)|l-2t| 

3 1/2  j  3l/2  j  jl/2 

has  maximum  value  at  t— j  +  (or  j-  -g— ),  hence  we  get  (52)  for  i— 0.  Similarly 

u,p-u np  x-x„  x.-x  „  .. 

T’(u(x);p)  -  ^UjP  — 5+  ^P  ^-}(»-P)/P 
on  h  h 

•  u’{  1  +  yyu(-j)”[(Xj-x)2+(x0'X)2-4(Xj-xXx-x0)J  }  +  0(hs) 
and 

^«P*unP«  x-xft  X.-X  tin  \ i 

T"(«(k);p)  -  (1-rtM-iftV  -2+ 

X.+X0  y>2  U  - 

- »”{ i  *  (“y-2-  -  + ^h1  + 

!^^^p*1u’)’[(x1-x)u1P+(x-x0)u0P)  +  0(h5)}. 

Hence 

T’(u(x);p)  -  u’(x)  -  T^y”[h2  -  e(x,-xXx-x0)]  +  0(h3), 


(54) 


(55) 
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T”(u(x);p)  -  u”(x)  -  •  (^4^-  ’  *frHT  +  0(h2).  (56) 

2  ®  u’ 

(55)  and  (56)  lead  to  the  bounds  (52)  for  i=l,2,  respectively.  In  order  to  prove  the  last 
estimate  of  (52),  nsing  (54)  we  find 

I|T-.||V  -  /0‘[T(tt(x) )  *  a(x)l2dx  -  ^A2^. 

Remark:  It  is  interesting  that  the  coefficients  C0,Cj,C2  in  (53)  are  just  as  same  as  Cj,C2,C3 
respectively  in  an  error  estimation  which  is  for  the  cubic  Hermite  interpolation  in  the  case  of  ufC4. 
e.g.  [23]  (theorem  2.21).  Hence,  it  is  reasonable  to  call  this  kind  of  piecewise  interpolation  a  quansi- 
cubic-hermitian. 

Corollary.  If  ueC3,  then  the  main  orders  in  (52)  still  keep,  however,  the  constants  before  the 
orders  are  need  to  change  now. 

We  may  extend  the  interpolation  form  (46)  further.  In  general,  suppose  {F^}  is  a  sequence  of 
piesewise  one-to-one  mappings  in  the  subintervai  [x^j,Xj]  respectively.  For  a  fixed  j,  say  j  «=  1,  let 
F—Fj,  we  define 

T(u;p)  -  F{tF‘*Uj  +  (l-tr1^)  (57) 

where  the  notation 

F-luj  -  F*1u(x)|  (i  -  0,  1),  t  - 

i  h 

It  means  that  now  the  piecewise  linear  interpolation  operator  does  operate  on  the  map  of  F*lu 
instead  of  on  u  itself  directly.  In  this  sense,  (46)  is  merely  an  example  of  (57)  where  Fu  “  u*/p. 

Since  T(u(t);F)  — ;  u(t)  at  two  ends  t— 0  and  t»l,  it  leads  to 

Lemma  14:  Let  u,  Ft  Cki'1  (  x  1  /  where  k"*$  or  4,  F1**  ie  any  one-to-one  mapping, 

eay  ^  >  0,  then 

•  i  h2  d2F'!u(x)  , 

lir>T(«)  •  ^  II— +  o(k"), 
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d  ,  ,  h  d2F*u(x)  .  . 


Theorem  15:  Suppose  the  hypotheses  of  Lemma  14  hold,  then  there  are  remainder 
formulas  and  error  estimates 

d2F(T  F£u(x)) 

u(x)  -  T(u;p)  —  ^x-xHXxj-xXu” - - - )\(, 

d2F'*u(x)  dF1n(x)  ,  dF',o(x)  , 

u(x)  -  T(u;p)  —  J(x-xHXxj-x){u”tf  )  +  u’2 — — ( — -j—  Y*\£ — -j— -  ft,,}, 


(  Xj.i  <  tv  <  X:  ) 

lnt»)  -  «lloo  Sj  Ih-^lloo^l.-F',  +  °(kk). 

s.||H^l|L^l!||_rHl  +  0(kk>). 


(58) 


(59) 


Proof:  The  first  remainder  is  obvious.  The  second  one  needs  differential  formulas  in 
implicit  form 

d^TF'Mx))  d  dp  .dTF'Mx)^  d2F-£n(x)  ,dF  »u(x)  ^dTF'Vxh 

— t~2 - -  a^ar^-r1  (  d»  r - tt~  (~ ar^  r  (  ix  * 


dx2  di'dt,,s-r*'  dx  '  dx2 

value  theorem 

dF*1  u(x) 


and  the  Mean  value  theorem 
F1u,-F1 


"  1  -  ~V  - 

Hence,  (58)  is  proved.  From  Lemma  14, 
h2  d2F£u(x) 


F*T(u)  F*u  +  — |h 


IL  +  o<kk) 


8  "  dxs 
Being  monotony  increase  of  F1 

h2  d2F‘*u(x)  . 

T(u)  F(Flu  +  ~.||__||oo  +  0(hk)} 

using  the  Taylor  expansion  completes  the  proof  of  the  theorem. 

F  urthermore 

Theorem  16:  Let  u,  F  t  C*  [zff  x}],  if  there  exists  a  (  t  (xQ  ,xt)  such  that 
d2F*1u(x)| 

dx2  lx“e  ~ 


0, 


h 

3 

»' 

F 

I 

s 

1 

b 

I 

»■ 

j 

4 

4 

< 
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then 

l3  d^'^x)  jp/_\ 

iit(u)  .  uiu  <l  V  + w4)-  m 

Besides,  if  (  —  (x^zJ/2,  the  coefficient  8 ’  in  the  dominotor  of  (60)  ean  be  improved  by 

•16’. 

Corollory.  T(u)  ■*  u  for  all  x  if  and  only  if 


dVMx) 


i.e.,  u(x)  —  F(C0-f-Cjx),  where  C0,Cj  —  constant. 

Observing  that  the  above  piecewise  interpolator  functions  (40)  and  (57)  all  only  belong  to  C°, 
nevertheless,  we  have  also  designed  a  piecewise  sublinear  positive  interpolator  function  which 
belongs  to  C1  [22]. 


5.  An  application  to  a  singular  perturbation  boundary  value  problem 
Consider  the  following  boundary  value  problem 
Lu  -  -<  u"  +  p(x)  u’  +  q(x)  u  -  f(x), 

u(0)  —  u(l)  —  0  (61) 

where  t  is  a  small  positive  parameter  and  p(x),q(x)  and  f(x)  are  so  sufficient  smooth  that  their 
derivatives  until  second  order  are  uniformly  bounded  for  all  x  in  [0,1]  sad  for  all  e  >  0,  besides,  p(x) 
p*  >  0  ,  q(x)  ^max(0,  p’(x) )  on  [0,1]. 

Let  Hm  be  Sobolev  space  of  m-order  with  the  norm 

IMI,  -  (  (D'o)2  )‘/! 

and  a(u,v)  be  the  unsymmetric  bilinear  form 

a(u,v)  — >  /0*  (  <uV  +  pu’v  +  quv  }  dx  (62) 

With  these  notations  the  weak  solution  of  (61  )  can  be  written  as  ;  Find  u  e  H°,[0,1]  so  that 

a(u,v)  ■»  (f,v)  for  all  v  «  H°j[0,l]  (63) 

where  H°,[0,1]  ~  {  v|  v  <  H,[0,1]  and  v(0)  -  v(l)  -  0  } 

Existence  and  uniqueness  of  solutions  to  (63)  follow  from  Theorem  5  using  the  following  Lemma: 
Lemma  17:  /II )  There  eziete  a  constant  C> 0  which  is  independent  of  t  ouch  that 

|a(u,v)|  C  Hull,  t  Hvllj  for  all  u,v  «  H°, 

|a(u,v)|  C  Hull,  t  IMIlit|1/(  for  all  u,v  e  H°,  (64) 

and 

N*»u)|  C*1  ||«||lf*,  for  all  u  t  H°,  (65) 

where 


Hull,.  -  (/„'  <«.'*  +  U=)d*>V2 
Hull,  .  ,/.  -  (/„'  (<«*  +  i«!)dx)'/* 


m 

(#7) 
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Now  we  apply  the  generalised  Galerkin  method  described  in  Section  3.  Because  the  singularity  of  the 
solution  u(x)  of  (61  )  is  only  near  x**l,  the  width  of  the  boundary  layer  in  which  u(x)  has  large 
derivatives  is  less  than  k  times  t,  where  k  is  a  constant  no  matter  how  t  is  small,  and  on  [0,1-ke]  u(x) 
and  its  some  first  derivatives  are  uniformly  bounded. 

Let  Ah  denote  a  partition  of  the  interval  |0,l]  into  N  subintervals  [xj _j»Xj],  j=l,2,...N  with 
x0=0,XjN-=l.  For  convenience,  we  will  consider  only  the  case  of  uniform  mesh  :  x^-x.  j  «■=  h,  j  *■= 
1,2, Associated  with  Ab  we  have  two  subsets  with  same  freedoms  of  H°j[0,l],  one  is  the  usual 
piecewise  linear  space  Pb  ,  another  is  called  SPjb  which  is  defined  by  that  if  u#b(x)  t  SPjb,  then  for 
Xj.!  ^.x^Xj,  t  =  (x  -  xH)/h, 

uj-i(1~t')  +  V  if  luj  ‘  uj-il/h  <  dl 

ush(x)  -  {  (68) 

(«j.1+c){(«j+c)/(“i.l+c)}t  -  *  Otherwise 

where  c  is  a  parameter  to  be  such  chosen  that  it  makes  the  formula  to  be  well  defined  and  to  get 
better  approximation  for  the  special  problem,  dl  is  a  controllable  constant. 


For  a  fixed  u(x),  the  interval  (0,1)  now  divides  into  two  subintervals  :  [0,1]  «*  If  +  Ig,  where  lr 
will  be  be  called  regular  on  which  the  first  derivative  of  u(x)  is  bounded  by  a  control  number,  Ig 
-  singular  subinterval  in  which  u’(x)  could  be  very  large. 

Being  Theorem  12,  for  fixed  c  and  dl,  SPtb  consisted  by  all  admissible  elements  of  (68)  is  a 
sublinear  set  of  H°j,  it  is  differs  from  the  corresponding  linear  space  Vb  only  where  the  element  has 
large  first  derivative. 


Let  (v.)  be  the  'roof'  basis  of  the  test  function  space 
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(x  "  xj-l)/h  xj-i^.x^-xj 

Vjh(x)  «=  {  ( j  —  1,2,...,N-1  )  (69) 

(Xj+l*x)/h  5Cj^.X^.Xj+, 

For  the  sake  of  simplify  we  first  suppose  that  the  coefficients  p  and  q  in  (61)  are  constant.  In  order 


to  get  the  integration  (62)  we  need  the  following  Lemma  which  can  be  convinced  by  part  integration 
Lemma  18:  For  ab  >  0, 


In  particular 

1  b  -  a 

I  _ - (b - }. 

1  Log(b/a)  Log(b/a) 

There  are  some  inequalities  in  [21]  about  IQ  and  1,  which  will  be  used  later: 
Lemma  19:  Suppose  a,b  >  0,  then 


(ab)1/2  ^.I0 

-(ab)1/2  minfLa'1/4^/4)  ^.1,  - max(l, 

2  4 

with  ’■*’  iff  a  *=  b. 


b+(ab)1/2 
b  +  a  ^ 


(70) 


(71) 


The  corresponding  integral  of  linear  interpolation  to  is 

a  +  (k+l)b 


LIfc  -  /01[a(l-t)+bt]tkdt 


Therefore  we  have  estimates 


(k+lXk+2) 


O^LI0-I0^^(b,/2-a,/2)2I 

j^(b1/2-a1/2Xb1/2-2a1/2)  ^  Llj  -I,  ^  J(b1/2-a1/2X2b1/2-a1/2)(  (b^a>0), 


*  yr^a-b)  ^  LI,  -  I,  J(a1/4-b1/4Xa3/4+a,/2b1/4+a1/4b1/2-  2b3/4),  (a£b>0), 

SUP  !LIo  *  i0l  ™  V  8UP  lLIl '  !ll  “  4- 

0<»,b£.l  w  w  *  0<»,b£.l  1  1  » 


Integrating  (62)  from  x.  ,  to  Xj  yields 


(72) 


/01(K+PtXc+«j.,)I*t(c+«j)tLo8 - J-dt  +  bq  /0,{(c+Oj.i),'t(c+«j)t  *  c}tdt 


-(Uj-U;  .)  +  p{c  +  Uj - -  }  + 

b  J  J‘  1  LogUc+up/te+Uj.,) ) 


bq{*  |  + 


Uj  +  C 


Vw 


2  Log((c+Uj)/(c+Uj_,)  )  (Logftc+Uj^/fe+Uj)))2 

Similarly,  integrating  (62)  from  Xj  to  Xj+1 

a(ugh,vj+h)  - 

J0'(-  j+pO-1)  Xc+uj),*t(c+uj+1)tLog - —  dt  +  hq  /0l{(c+u.),*t(e+u.+1)t  -  c}tdt 


ui+i*ui 


r<uruj+i)  *  p^c  +  uj)  •  7~~7( — l~T7, — ;  J  + 
h  1  1  1  Log((e+uj+1)/(c+Uj) ) 


U;  +  C 


hq.<-  £  + - “J—1 - - 'lZi±l -  j. 

1  2  Log((c+Uj)/(c+nj+1) )  (Log((c+Uj+1)/(c+Uj)))2 

For  (Xj  j,  Xj]  (  If,  a  straightforward  computation  yields 

a(u.h-vjh)  “  s12Vh*Vi1  +  4p(ui+ruj-i)  +  |favi+4uj+uH) 

and  for  [Xj  , ,  Xj]  <  Ig 

a(ugh,Vjh)  —  a(ugh,Vjh)  +  a(ugh,vj+h) 

-  t(2u  -u  r«.+1)  +  p  { - - - -  }  + 

F  J  J*  J+  I^o8((c+U; , , )/(c+u;) )  Log((c+u.)/(c+u. ,)  ) 


gUc+u:)/(c+Uj . )  Log((c+u.)/(c+u. ,)  ) 


(Uglic+Uj^jj/tc+Uj)))2  (Log((c+Uj_l)/(e+uj)))2 


I 
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a(ugh,vjh)  -  j(2uj-uj.1-uj+1)  +  5P(uj+1-tij.1)  +  5^uj+,+4uj+uj-,)  +  Sj 
where  gj  is  the  difference  of  the  right  parts  between  (73)  and  (74). 

Denote  ^  ■»  a,  substituting  (73  )  and  (75  )  into  the  generalised  GaJerkin  method,  i.e. 
a(ugh,  v.h)  —  (ftv>)  for  j  -*  1,2,...N-1 

leads  to 

(f,Vjh)  if  j  <  Ir 


*  { 


(f,^6)  •  g(U“H,  Uhj(  Uhj+1)  ifjcl, 
where  the  left  side 

LhUh  _  .(„+  1.  \)Ufc-,  +  (2o  +  -  (o.  I  ■  ^)UkH 

which  is  exactly  the  same  to  the  scheme  from  usual  piecewise  linear  subspace. 

With  matrix  form  it  can  be  written  as  the  special  form  as  (33). 

A  U  -  d  +  Q(U) 

where  A  is  a  tridiagonal  matrix  A  »  (  a. .  ) 


-(°  +  f-£q)  *>i 

2h 


(75) 


(76) 


(77) 


(78) 


(79) 


oiJ-(2o+J!iq  i-j 
•(a-f'jq)  i<j 

Denote  the  determinants  of  the  first  j  and  the  last  N-i  principal  determinants  of  A  by  Dj  and  D;  N  l 
respectively,  set 

t  -  Esj,  f.  „ ,  -  Pi+W. 

®  n  d 

n  Uj,N-l 

Due  to  the  recursion  formula 

(2o  +  |^q-(a- J-gqXo  +  f-gq)^}'1. 
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[ 

i 


» ■ 

therefore 

c 

Lemma  20:  If  a  —  ^  ^  ^  then 

*n  fHo+f-gq}'1,  for  all  n  <  N-l. 

0nN.j  £.{<*  +  \  -  gq}'1.  for  all  n  <  N-l. 

—  (80) 

1 

Meanwhile,  we  have 

Theorem  21:  When 

°  +  S*1 

(81) 

1 

A'1  —  (  a '■*.  .  )  is  a  good  discrete  Green  function  in  the  following  sense:  A'1 
hj 

is  non- 

■ 

negative  and 

iriii  or  i°'Vi  ifi<i 

(82) 

'« 

Proof:  In  fact,  in  this  case  A'1  «■*  (o'1- ;) 

*tJ 

o-'u  -(<«  +  !■  kiDHDN.l,/DN.,  « 

(«-f-jOi-iD.1DN.1.i/DN.1  if*J 

(83) 

Since  A-1  exists,  (78)  can  be  written  as 

» 

U  -  A  *(d  +  Q(U)  ) 

(84) 

Now  we  look  for  an  estimate  of  ||A*IJ(Q(U)  )||,  where  J(Q)  is  the  Jacobi  matrix  of  Q.  The  main  idea 

of  the  derivation  is  the  same  to  oar  another  paper  [21]  in  which  the  scheme  based  on  a  second  order 

semi-linear  numerical  differentiation  formulas  has  the  same  form  (84)  with  slight  different  A  and 

Q.  Thus  we  only  need  to  explain  the  outline  of  proofs  which  are  different  here.  First,  we  prove  that 

the  following  important  ’semi-linearity’  of  Q  defined  in  (21): 

Lemma  22:  For  Q(u)  defined  by  the  difference  between  the  linear  scheme  (73)  and  the 
semilinear  scheme  (74),  there  exits  the  following  identities: 

(J(Q(u)  X«+c)}j  -  {Q(u)}jf  if  j  <  N-l.  (85) 

Proof:  Denote  the  nonlinear  term  of  third  term  in  the  right  of  (74)  by 

F("-1'"«"1)  - e  +  “o  •  ,c+“o)la,i(i;+l«tiAc+;.l)) +  Mre'+g/it+u,))1  + 


I 
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“l-^  ,  ”o  ~  tt-l 

(Log((c+u1)/(c+u0)))2  (Log((c+u0)/(c+u  ,)))2 

3F  _  .(c+u  —  *  _ l- _ 

C+  °(Log((c+Uo)/(c+u  j)))2  (Lof((c+u0)/(c+u  j)))2 

2(V”|)/(C+U.l) 

(LogHc+u^c+u  t)))3 


3F  f- 

1  t 


1  1 

-  + - - —  ]  + 

ogMc+u^/fc+u^)  Log^e+u^/fe+u  j)) 


(c+uo)(: 


l/(c  +  1^) 


l/(«  +  «o) 


(LogUc+u^/fc+u  ,)))2  (Logftc+u^/fc+u,)))2 


(LogUe+u^Ae+u  ,)))2  (Log((c-l-a0)/(c-t-n1)))2 

.  i_  l  v*  +  va-»  , 

c+u0  (Log^c+u^/tc+u,,)))2  (Logtfc+U^Ac+U  j)))2 

dF  _  (c+  )  l/(c  +  °i)  (  1 _ 

^  +U°  (Logftc+u^/fc+u,)))2  (LogHc+u^/fc+u,)))2 

2(”f”0)/(c+tti) 

(LogUc+u^/lc+u,)))3 

Hence 

5^(c+a-i)  +  5^c+uo)  +  5^c+ui)  -  F(»  1-Vi)' 

It  has  been  proved  in  (21]  that  the  second  term  in  the  right  of  (74)  satisfies  (85),  due  to  the 
linearity  of  the  ’semi-linear'  relation  we  get  (85). 

Since  the  singularity  is  only  near  x«l  and  the  width  of  the  boundary  layer  is  less  than  let,  using  the 
inequalities  (72,  (82),  (81)  and  (83),  similarly  to  (21]  a  straightforward  computation  yields 
J(Q(u)  Xu+c)  -  (0 . 0,Qn,...,QN-2,Q*N_|}, 
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{A-1J(Q(u)  )u};  —  (Tj  +  hqr; 


where 


*  |£<  (®  +  f  -  <  ?>  -  s**  <  (®  +  5  ■  s^>b<iri  < 


Therefore  — 

Theorem  23:  When  (he  mesh  size  condition  (81),  i.e., 

k&|^J+lj+§({?0W-,’H»-!<{?->!>  (»» 

holds  as  well  as 
3p 

(87) 

8q 

then  the  mapping  AmlQ(u)  is  contractive,  in  the  meantime  the  semi-linear  system  (78)  can 
be  solved  by  the  following  convergent  'simple  ’  iteration 

A  U<°>  —  d 

A  UW  -  d  +  Q(U<kl>)  (k  -  1,2,...)  (88) 

2f 

Remark:  When  e  is  small, in  practice,  the  mesh  condition  (86)  can  be  simplified  by  h  < 

Now  we  consider  error  estimations.  Let  u  be  the  true  solution  of  (61),  there  exist  a 
decomposition  [11] 

u(x)  -  7{W(x)  +  Z(x)} 

W(x)  -  -  x  -  (1-xKPt1)/*  (89) 

where  <7  —  lim  lim  <u’(x)/p(l)  is  a  constant  bounded  uniformly  for  all  0<(<1,  and 
«-  >  0  »-  >  1 

|Z(x)|£.C,  |Z’(x)|£.C,  |Z”(x)|^.C{l+I^,'*)/f} 

C  is  a  constant  independent  of  <,  and 


Set  c  in  (89)  equal  to  7  which  can  be  found  in  computing  test.  We  proved  in  [21]  that 
Lemma  24:  Let  v  be  the  (rue  solution  of  (61),  if  h  and  t  are  of  the  same  order,  then 

ll«(j,IL  -  O(h-i),  (j  -  1,2,...) 
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H«"  •  3„ll„  -  W).  IK"”  -  £.>11.  “  0(h'2'- 

where  e  *  lim  lim  tu'(x)/p(J). 


(90) 


Being  (49)  and  (50),  hence,  for  the  interpolation  function  uhj  of  u(x)  in  SPjh  we  have  error 
bounds 

ll»h,  -  «IU  -  0(k).  II»V  “'ll»  -  °<1>'  (»1) 

Moreover,  since  the  width  of  the  boundary  layer  is  the  same  order  of  t,  keeping  h  as  the  same  order 
of  t  too,  it  leads  to 

ll»h,  -  "II,  -  0(b3/=>,  Ha",  -  all,  -  CKh'/2), 

ll"kj  -  "II, ,<  —  0(b).  Il"11,  -  "Hi.,,,/,  —  0(b).  (92) 

Let  Hj  and  H2  be  the  Hilbert  space  with  the  norm  (00)  and  (07),  respectively.  Using  Lemma 
17  and  (31)  yields 

II*  -  *%,,  ^  ^ II*  *  uJhlll,e,l/<  +  IK*  -  «bH®  *  Alll,e,l/«  +  IK®  -  ®JhM»  *  «Jh)llli,(,i/f} 

2  h 

where  the  subscript  I  denotes  the  interpolation  in  the  test  space  Vn  -  piecewise  linear  function 

subspace.  On  the  right  side  of  the  above  inequality,  the  first  term  is  the  major  one,  others  are  of 

higher  power  of  h.  Hence,  being  (92),  we  get  the  main  error  estimation  for  the  scheme  (70) 

Theorem  25:  If  the  mesh  size  condition  (86)  holds,  then 

ll».h  -  "111.,  —  0(h>.  (93) 

where  coefficients  before  powers  of  h  are  uniformly  bounded  for  ail  small  c  satisfying  (87). 


Applying  the  Taylor  expansion  and  using  the  equation  (01)  itself  and  (72),  substituting  the  true 

solution  u  into  the  scheme  (77)  yields 

h3 


LhUj  —  (f,Vjh)  +  0(h2)  +  +  +  9®”}  +  - 

if  j  t  lr,  and 
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L^o.  —  (f,v>)  +  0(h2)  -  gjtuj.pUj.u^j)  +  Trj(u), 


B 


l 


i 

i 


.>2 


Tr-  -  -{P«(3)  -  p£fu  )’  +  2qu”)  +  ... 

12 

Using  (00)  and  noting  the  fact  that  the  width  of  the  boundary  layer  in  only  ke,  similar  to  [21]  we  get 
l|A-,Tr(x)||eo  -  0(h)  -  (94) 

furthermore 

!Kfc  -  “IL  -  °0>). 

Similarly 

ll*,k  -  »ll0  -  0(hi  s),  lln.*1  •  nil,  -  0(k«). 

Summarizing  the  above  results,  finally  we  obtain  the  following  theorem  of  error  estimations 
Theorem  20:  For  email  e  satisfying  (81),  when  the  mesh  size  condition  (86)  holds,  then 
the  generalized  Galerkin  method  on  the  subset  (16)  has  one  more  order  of  precision  than 
its  corresponding  scheme  of  piecewise  linear  sutspaee,  i.c.,  in  this  case  we  have  there  exist 
constants  and  C’^  which  are  uniformly  bounded  for  all  small  e  such  that 

ll»,h  -  "llo  S.C0h15,  Hu,"  -  n||, 

llV  -  «ll„  £  Coo11'  IKh  •  "’llo,  £  C’„-  (05) 

In  the  case  of  the  general  variable  coefficients  p  and  q,  it  can  be  proved  that  the  above 

conclusion  still  holds  true  for  small  e  if  two  extra  requires  are  satisfied: 

2  1  1 

h  <  Iip||  *’  **d  ^qj*,+4qi+qj+,)  j  “  1A*“  (96) 

the  last  one  is  a  discrete  form  for  the  elliptic  condition  of  afxl>p,fx). 


As  a  matter  of  fact,  we  only  need  to  point  that,  being  smoothness  of  p  and  q,  substituting  their 
piecewise  linear  interpolations  into  the  integral  form  (62),  (73)  becomes 

•<«>>)  -  H^j-Uj-rV,) +  ^uj+i<2Pj+Pj+i)+nj(Pj-rPj+i)-Bj-i(2Pj+Pj-i)] + 

^«j+i(qj+qj+i)+«j(qj.,+6qj+qj+,)+,ij.,(qj+qj.,)]  +  o(h2) 


(97) 
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1 


I 


( 


The  related  tridiagonal  matrix  A  ■»  (a. .)  in  (70)  now  is 
*(«  +  g(2Pj+Pj-i)  *  T2^j+qj-i)l» 

%  “  <  20  -  ^Pj+rPj-i) +  i-j  (®8) 

-[o  *  ^Pj+Pj^.,)  -  j^qj+qj+i)].  ‘-H  - 

The  rest  derivation  is  similar  to  (21],  we  omit  it  in  detail. 

For  higher  order  schemes  based  on  the  interpolation  described  in  section  4,  the  similar  analysis 
can  be  also  done. 
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6.  Numerical  Results 

Id  this  section,  we  give  two  examples  to  show  how  well  the  numerical  results  match  the 
conclusions  in  Theorem  25  and  26.  In  the  following  Tables,  the  notation  N  ■*  1/h,  SL  and  L 
represent  the  subset  scheme  (76)  and  its  corresponding  linear  scheme,  respectively,  Er(Max)  is  the 
maximum  error  with  sign  of  the  discrete  solution  and  it  has  occupied  on  the  node  xM,  Er(H  l,eps), 
Er(H0)  and  ErfHl)  are  the  approximation  values  of  errors  in  Hj  (,H0  and  Hj,  respectively.  CPU  -  the 
CPU  time  in  terms  of  seconds.  The  Fortran  program  was  run  in  double-precision,  on  a  DEC-System 
2060  computer.  The  iterative  error  for  (88)  is  equal  to  10'5. 


EXAMPLE  1.  A  linear  singular  perturbation  problem  with  constant  coefficients 
Lu  —  -  eu”  +  u’  +  (l+«  )u  *■  f(x),  in  (  0,1  ) 


u(0)  —  u(l)  —  0 

where  f(x)  «■  (l+<  X*-b)x  -  «a  -  b,  a  —  1  +  e'^1+<  ^*,  b  «  1  +  e*1,  with  true  solution  (see.  Figure  1) 
u(x)  ■*  e^1+f  +  e"*  -  a  +  (a-b)x 

In  our  case,  set  the  constant  c  —  1  in  the  scheme  (68),  see  [21]. 

The  results  listed  in  Table  1-4  (or  Figure  2-4)  show  that : 

1.  The  iteration  of  (88)  monotony  converges  if  the  ratio  h/e  <  2  ,  the  results  match  with 
the  theoretical  analysis  above,  the  SL-  scheme  is  much  better  than  Lrscheme  with  little 
more  CPU  time  cost  (  about  20%  for  small  c  )  for  the  same  mesh  size  h. 

2.  When  2  ^  h/c  ^  2.25,  the  iteration  seems  still  convergent,  but  oscillation  is  occupied 
now,  and  the  error  is  getting  more  than  the  above  estimates,  CPU  time  is  more,  too. 

3.  If  the  ratio  increases  again,  the  iteration  (88)  does  not  converge. 

4.  For  a  given  level  of  accuracy,  the  CPU  time  costs  much  less  using  the  SL-scheme  than 
using  the  L-scheme,  and  more  small  <  there  is,  more  advantage  the  SL-scheme  has.  For 
instance,  given  an  admissible  maximum  error  at  knots  0.005,  their  CPU  time  ratio 
are  about  0.3  :  1.1  and  3  :  15,  for  e  ■>  0.01  and  0.001,  respectively. 


EXAMPLE  2.  A  semi-linear  singular  perturbation  problem 
Lu  -  tu”  ■+  p(x)u’  +  q(x)u  ■■  f(x,u),  in  (  0,1  ) 


f(x,u)  -■  a  -  b  -  (l+<  ){e*x  •  u  + 


u+a-fa-bjx-e' 
a~l+e^1+<  )/*,  b—l+e'^c  ■»  e2^1+t  H1**)/*, 

p(x)«  l,q(x)  —  l+< 

with  the  same  solution  as  example  1. 


In  the  semi-linear  case,  the  advantage  of  SL  scheme  over  L-scheme  is  more  obvious  than  in 
linear  case,  the  results  of  S L-scheme  still  match  the  Theorem  25  and  they  are  much  better  than  L- 
scheme  with  same  conditions  to  obtain  higher  accuracy  and  save  computer  time  both  (see  Table  5-7, 
or  Figure.  5-7). 
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Table  1-1 


SL:  h/e  •  1.5 


N 

xM  ) 

Er(Nax) 

Er(Hl.eps) 

Er(BO) 

Er(Hl)  I 

1 

CPU 

25 

0.920  1 

-0.5837D-02 

0.35080-01 

0.1740D-02 

0 .21460+00 1 

0.09 

50 

0.960  I 

-0.6023D-02 

0.2079D-01 

0.1239D-02 

0.17980+001 

0.16 

100 

0.970  1 

-0.4562D-02 

0.1555D-01 

0.7205D-03 

0.19020+001 

0.50 

200 

0.985  1 

-0.1810D-02 

0.8410D-02 

0.2019D-03 

0.14560+001 

1.07 

400 

0.993  1 

-0.1239D-02 

0.5495D-02 

0.10310-03 

0.1346D+00! 

2.17 

800 

0.996  I 

-0.5259D-03 

0.3000D-02 

0.32230-04 

0.10390+001 

4.47 

1600 

0.998  I 

-0.3363D-03 

0.1823D-02 

0.14830-04 

0.89280-01 1 

8.56 

Table  1-2  L:  h/e  -  1.5 


N 

xM  | 

Er(Max) 

Er(Hl,eps) 

Er(HO) 

Er(Hl)  1 

1 

CPU 

25 

0.960  1 

-0.8199D-01 

0.1216D+00 

0.1742D-01 

0.7371D+O0I 

0.04 

50 

0.980  1 

-0.8112D-01 

0.1183D+00 

0.1223D-01 

0.1019D+01 1 

0.06 

100 

0.990  1 

-0.8070D-01 

0.1167D+00 

0.8620D-02 

0.1425D+01I 

0.33 

200 

0.995  1 

-0.8048D-01 

0.1159D+00 

0.6084D-02 

0.2004D+01 I 

0.75 

400 

0.998  1 

-0.8038D-01 

0.1155D+00 

0.4299D-02 

0.282 7D+01I 

1.53 

800 

0.999  1 

-0.8033D-01 

0.1153D+00 

0.3038D-02 

0.3992D+01 i 

3.23 

1600 

0.999  1 

-0.8030D-01 

0.1152D+00 

0.2148D-02 

0.5641D+01I 

6.54 

Table  2  SL: 

h/e  -  1.75 

N  I 

xM  1 

Er(Max) 

Er(Bl.epa) 

Er(BO) 

Er(Bl)  1 

1 

CPU 

25  1 

0.920  1 

-0.1899D-01 

0.4066D-01 

0.4869D-02 

0.2671D+O0I 

0.11 

50  1 

0.960  1 

-0.6682D-02 

0.2341D-01 

0.1254D-02 

0.2187D+00I 

0.25 

100  1 

0.980  t 

-0.4199D-02 

0.1643D-01 

0.6130D-03 

0.2171D+00I 

0.51 

200  1 

0.990  I 

-0.1913D-02 

0.1043D-01 

0.2120D— 03 

0.1950D+00I 

1.11 

400  I 

0.993  1 

-0.1162D-02 

0.6063D-02 

0.8863D-04 

0.1604D+00! 

2.48 

800  I 

0.996  1 

-0.7638D-03 

0.3817D-02 

0.4342D-04 

0.1428D+00I 

4.93 

1600  1 

0.998  I 

-0.3355D-03 

0.2080D-02 

0.1374D-04 

O.llOlD+OOl 

10.03 
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Table  3-1 


SL:  h/e  “  2.0 


N  1 

xM  i 

Er(Max) 

Er(Bl,eps) 

Er(BO) 

Er(Bl)  1 

1 

CPC 

25  1 

0.920  | 

0.6469D-02 

0.2969D-01 

0.1423D-02 

0.2097D+00I 

0.18 

50  | 

0.960  | 

-0.7583D-02 

0.2592D— 01 

0.1366D-02 

0.25880+001 

0.36 

100  1 

0.970  I 

0.3605D-02 

0.1359D-01 

0.3779D-03 

0.19210+001 

0.66 

200  1 

0.990  | 

-0.1235D-02 

0.1032D-01 

0. 10330-03 

0. 20640+00 1 

1.47 

400  I 

0.990  | 

0.1452D— 02 

0.5153D-02 

0.7843D-04 

0.14570+001 

2.93 

800  t 

0.996  I 

-0.6512D-03 

0.4100D-02 

0.3331D-04 

0.16400+001 

5.81 

1600  I 

0.997  I 

0.8136D-03 

0.1324D-02 

0.3186D-04 

0.74850-01 I 

12.21 

Table  3-2  L: 

h/e  -  2.0 

N  1 

xM  i 

Er(Max) 

Er(Bl,eps) 

Er(HO) 

Er(Hl)  1 

. . .  1 

CPU 

25  1 

0.960  | 

-0.13660+00 

0.16470+00 

0.27530-01 

0.11490+011 

0.03 

50  I 

0.980  | 

-0.13600+00 

0.16140+00 

0.1939D-01 

0.16020+011 

0.21 

100  1 

0.990  | 

-0.13560+00 

0.15970+00 

0.13690-01 

0.22500+011 

0.37 

200  I 

0.995  I 

-0.13550+00 

0.15890+00 

0.96680-02 

0.31710+011 

0.77 

400  1 

0.998  | 

-0.13540+00 

0.15840+00 

0.68330-02 

0.44770+011 

1.55 

800  t 

0.999  | 

-0.13  540+00 

0.15820+00 

0.48300-02 

0.63260+011 

3.35 

1600  I 

0.999  I 

-0.13 54D+00 

0.15810+00 

0.3415D-02 

0.89420+011 

6.61 

Table  4  SL 

h/e  -  2.25 

N  1 

xM  | 

Er  (Max) 

Er(Hl,epe) 

Er(BO) 

Er(Bl)  1 

1 

m 

25  1 

0.920  | 

0.1618D— 01 

0.2526D-01 

0.3461D-02 

0.1877D+O0I 

0.17 

50  I 

0.940  | 

0.2635D-01 

0.2103D-01 

0.5016D-02 

0.2167D+00I 

0.48 

100  1 

0.980  I 

-0.4804D-02 

0.1985D-01 

0.5968D-03 

0.2976D+00I 

0.86 

200  I 

0.990  I 

-0.1223D-02 

0.1143D— 01 

0.9819D-04 

0.2424D+00I 

1.46 

400  j 

0.985  | 

0.1082D-01 

0.7342D-02 

0.9483D-03 

0.2184D+00I 

3.86 

800  | 

0.988  I 

0.2192D-02 

0 .42710-02 

0.9035D-04 

0.1812D+00I 

8.98 

1600  | 

0.997  | 

0.4134D-02 

0.4955D-02 

0.2310D-03 

0.2970D+00I 

16.28 

1 
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Table  5-1 


SL:  h/e  -  1.5 


N 

T 

xM 

1 

Er(Max) 

ErOU.ept) 

Er(EO) 

Er(Hl)  1 

1 

CPU 

25 

1 

0.920 

1 

-0.1930D-01 

0.3845D-01 

0.4679D-02 

0.23370+001 

0.34 

50 

1 

0.960 

1 

-0.5220D-02 

0.1843D-01 

0.104004)2 

0.15930+00 1 

0.98 

100 

1 

0.970 

1 

-0.4529D-02 

0.1436D-01 

0.6564D-03 

0. 17570+00 1 

2.42 

200 

1 

0.985 

1 

-0.1742D-02 

0.7681D-02 

0.17870-03 

0.13300+001 

5.40 

400 

1 

0.993 

1 

-0.1201D-02 

0.5141D-02 

0.95040-04 

0.12590+00 1 

11.58 

800 

1 

0.996 

1 

-0.5050D-03 

0.2820D-02 

0.2971D-04 

0.9767D-01 | 

23.48 

1600 

1 

0.998 

1 

-0.3337D-03 

0.1738D-02 

0.13970-04 

0.8515D-01 I 

47.13 

l. 


Table  5-2  L:  h/e  -  1.5 


N 

1 

xM 

1 

Er(Max) 

Er(Hl.epi) 

Er(BO) 

Er(Hl)  1 

-1 

CPU 

25 

0.960 

-0.8244D-01 

0.12170+00 

0.17730-01 

0.73770+001 

2.64 

50 

1 

0.980 

1 

-0.8125D-01 

0.11840+00 

0.12310-01 

0.10200+011 

9.63 

100 

f 

0.990 

1 

-0.8100D-01 

0.11670+00 

0.86850-02 

0.14260+01)  25.26 

200 

1 

0.995 

1 

-0.8059D-01 

0.11590+00 

0.610204)2 

0.20050+011 

56.46 

400 

1 

0.998 

1 

-0.8043D-01 

0.11550+00 

0.43050-02 

0.28270+011 

11.24 

800 

1 

0.999 

1 

-0.8035D-01 

0.11530+00 

0.30410-02 

0.39920+011 

23.49 

1600 

1 

0.999 

1 

-0.80310-01 

0.11520+00 

0.2149002 

0.56410+011 

48.78 

Table  6-1  SL 

b/e  -  1.75 

N 

1 

xM  1 

Er(Max) 

£7(81, ep») 

Er(HO) 

Er(Bl)  1 

1 

CPU 

25 

0.920  j 

-0.19680-01 

0.390104)1 

0.47240-02 

0.25620+001 

2.81 

50 

1 

0.960  | 

-0.601704)2 

0.20800-01 

0.10410-02 

0.19430+001 

1.14 

100 

1 

0.980  1 

-0.38790-02 

0.150304)1 

0.549404)3 

0. 19860+00 i 

3.26 

200 

1 

0.985  I 

-0.343104)2 

0.113604)1 

0.35000-03 

0.21240+001 

7.61 

400 

1 

0.993  I 

-0.114304)2 

0.566204)2 

0.808004)4 

0.14980+001 

14.67 

800 

1 

0.996  1 

-0.753404)3 

0.36250-02 

0.40760-04 

0.13560+001 

29.78 

1600 

1 

0.998  1 

-0.324604)3 

0.198804)2 

0.12900-04 

0.10520+001 

65.90 

Table  6-2 
h/e  ■  1.75 


N 

xX 

Er(Max) 

25 

0.960 

T 

-0.1163D+00 

50 

0.980 

1 

-0.1077D+00 

100 

0.990 

1 

-0.1104IH-00 

200 

0.995 

1 

-0.1076IH00 

400 

0.998 

1 

-0.10751H00 

800 

0.999 

1 

-0.1073D+00 

1600 

0.999 

1 

-0.1072D+00 

Er(Bl,eps) 


0.1453D+00 
0.1411D+00 
0.13 97D+00 
0.1387D+00 
0.1382D+00 
0.1380D+00 
0.1379D+00 


Er(BO) 


0.2433D-01 
0.1569IV-01 
0.1145D-01 
0.7839D-02 
0.5538D-02 
0.3908D-02 
0.2760D— 02 


Er(Hl) 


0.9480D+00 

0.1312D+01 

0.1842D+01 

0.2590D+01 

0.3655D+01 

0.5163D+01 

0.7297D+01 


2.73 

10.29 

26.24 

56.17 

119.92 

140.12 

496.64 


25 

50 

100 

200 

400 

800 


.920 

-0.1928D-01 

0.3572D-01 

0.4404D-02 

.960 

-0.7551D-02 

0.2358D-01 

0.123 8D-02 

.980 

-0.3771D-02 

0.1592D-01 

0.4465D-03 

.985 

-0.2508D-02 

0.8006D-02 

0.2353D-03 

.993 

-0.1181D-02 

0.6546D-02 

0.8375D-04 

.996 

-0.7179D-03 

0.4017D-02 

0.3538D-04 

.2507D+00  3 .42 

.2354IH00  10.95 

.2251D+00  25.85 
.1601D+00  55.86 
.1851D+00  116.21 
.1607D+00  236.06 


Table  7-2 


N  I  xX 

1 

Er(Kax) 

Er(El.epe) 

Er(HO) 

Er(Bl) 

CPU 

25 

0.960 

-0.1421D+00 

0.1653D+00 

0.2877D-01 

0.1151IH01 

3.28 

50 

1 

0.980 

-0.1337D+00 

0.1612D+00 

0.1907D-01 

0.1601IH-01 

11.33 

100 

1 

0.990 

-0.1255D+00 

0.1590D+00 

0.1264D-01 

0.2242D+01 

27.03 

200 

1 

0.995 

-O.13691H00 

0.1590D+00 

0.9775D-02 

0.3173D+01 

58.68 

400 

1 

0.998 

-0.1337D*00 

0.1583D+00 

0.6743D-02 

0.4474D+01 

124.43 

800 

1 

0.999 

-0.1323D+00 

0.1580D+00 

0.4715D-02 

0.6318IH01 

248.85 

Figure  1 


TRUE  S0LUTI0N  0F  EXAMPLES 


Figure  2  EXAMPLE  1:  -  -  1.5 
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C0MPARS0N  EXAMPLE  1,  H/EPS  =1.5 


1.0E+01 


l.OE+OO 


l.OE-Ol 


1.0E-02 


LINEAR  SCHEME  MAXIMUM  ERROR 
LINEAR  SCHEME  HO  NORM  ERROR 
LINEAR  SCtCME  HI  NORM  ERROR 
LINEAR  SCHEME  CPU  TIME  SEC. 
SEMIL INEAR  SCHEME  MAXIMUM  ERROR 
SEMILINEAR  SCHEME  HO  NORM  ERROR 
SEMIL INEAR  SCHEME  HI  NORM  ERROR 
SEMILINEAR  SCHEME  CPU  TIME  SEC. 


.i-UJillll 


Fifore  3 


C0MPARIS0N.  EXAMPLE  1.  EPS  =  0.( 


1  LINEAR  SCHEME  MAXIMUM  ERR0R  P0INT 

2  LINEAR  SCHEME  MAXIMUM  ERR0R 

3  LINEAR  SCHEME  FIRST  DERIVATIVE  ERR0R 

4  LINEAR  SCHEME  CPU  TIME  SEC. 

5  SEMILINERR  SCHEME  MAXIMUM  ERR0R  P0INT 

6  SEMIL INEAR  SCHEME  MAXIMUM  ERR0R 

ft  SEMILINEAR  SCHEME  FIRST  DERIVATIVE  ERR0R 
6  SEMILINEAR  SCHEME  CPU  TIME  SEC. 


Figure  4 


SEMI -LI NEAR  EXAMPLE  1.  H/EPS 


1  P0SITI0N  XM  0F  MAXIMUM  ERR0R 

2  MAXIMUM  ERR0R 

3  ER(H1»EPS) 

4  ER(HO) 

B  ER(H1) 

6  CPU  SEC0NDS 


Flgoxe  5 


C0MPRRIS0N  EXRMPLE  2.  H/EPS 


1 . 0E-Q5 


1  LINEAR  SCHEME  MflXIMUM  ERROR 

2  LINEAR  SCHEME  HO  N0RM  ERROR 

3  LINEAR  SCHEME  HI  NORM  ERROR 

4  LINEAR  SCHEME  CPU  TIME  SEC. 

5  SEMIL INEAR  SCHEME  MAXIMUM  ERROR 

6  SEMILINEAR  SCHEME  HO  NORM  ERROR 

7  SEMILINEAR  SCHEME  HI  NORM  ERROR 

8  SEMILINEAR  SCHEME  CPU  TIME  SEC. 


1.DE+Q2 


1.0E+01 


Figure  6 


C0MPRRIS0N.  EXRMPLE  2.  EPS 


1 . OE+OO ! 


1 . 0E-03 ' 


1.0E-04  i 


1 . 0E-05 


1.0E-06 

500. 


875. 


1250. 


1625. 


LINEAR  SCHEME  MRXIMUM  ERROR  POINT 
LINEAR  SCHEME  MAXIMUM  ERROR 
LINEAR  SCHEME  FIRST  DERIVATIVE  ERROR 
LINEAR  SCHEME  CPU  TIME  SEC. 

SEMIL INEAR  SCHEME  MAXIMUM  ERROR  POINT 
SEMIL INEAR  SCHEME  MAXIMUM  ERROR 
SEMI LINEAR  SCHEME  FIRST  DERIVATIVE  ERROR 
SEMILINE AR  SCHEME  CPU  TIME  SEC. 


Fignr*  7 


CPU  C0MPRRIS0N  EXRMPLE  2 

i - 1 - 1 — i — i — i — i — i — f 


1/H 

1  LINEAR  SCHEME  H/EPS  =1.5 

2  SEMILINEAR  SCHEME  H/EPS  =1.5 

3  LINEAR  SCHEME  H/EPS  =1.75 

4  SEMI LINEAR  SCHEME  H/EPS  =1.75 

5  LINEAR  SCHEME  H/EPS  =  2 

6  SEMI LINEAR  SCHEME  H/EPS  =  ? 


